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However, sandy deserts or sandy playas, with virtually no vegetation, might be an almost unique 
terrain type in this respect, that they reasonably approximate a Lambert surface when observed 
from the zenith. Vegetated areas or urban areas can hardly be regarded as planar, because most 
plants and buildings form sharp protrusions from the plane of the soil. Whether the protrusions are 
on a scale of a blade of grass or of a skyscraper is immaterial; in either case, the protrusions invalidate 
the Lambert Law for the surface reflectivity 4,5 . Thus, the analysis of 9 Landsat passes, all within 
one year, over a presumably unchanging surface (an ungrazed rangeland in Utah), showed variations 
in the reflectivity by about 25% in all the spectral bands, from the January and December lows to 
the June highs. The variations were largest in the 0 8-1 . 1 jam band (~30%) and smallest in the 
0. 5-0.6 Mm band. The effect can be explained as a systematic dependence on the solar zenith angle 
of the reflectivity to the zenith of the surface. The plants in this arid steppe form 10 to 20 cm high 
tussocks or bushes protruding from the soil 5 . Their shadows, which depend on the solar elevation, 
quite obviously reduce the observed reflectivity in winter as compared to summer, for observations 
at a fixed local time. (Landsat is in a sun-synchronous orbit, with the subsatellite local time at about 
09:30 h). The infrared to red reflectivities ratio of this site remains essentially constant throughout 
the year. The measured reflectivities of a playa, a few kilometers away, also remain approximately 
constant. 

The author's own interest in the problem started with the observations of a fenced-off area in 
the Sinai, where a rapid recovery of the natural vegetation took place after the entry of the grazing 
herds was stopped by the fence 3 . The plants consist of mostly grey or dark brown twigs. The green 
components are quite inconspicuous, except in the desert-bloom period following a rain. Each plant 
forms an isolated clump, with large interstices of bare soil between a plant and its neighbors, see 
Fig. 1, in order to draw moisture from a large area. In this study, we analyze the atmospheric effects 
that can complicate the task of monitoring such a desert-fringe area from a nadir-viewing satellite. 

Monitoring a non-Lambertian surface by radiometric imaging from a satellite can be rather 
difficult, inasmuch as there is no simple correspondence between the differences in the value of 
the reflectivity measured at different satellite passes, and the possible changes in the surface. 


than those over a Lambert plane for the same surroundings-to-object-pixel contrast and atmos- 
pheric conditions. However, the adjacency effects are highly variable, because the effective con 
trast for our plane w th dark protrusions is a function of not only the surface parameters but 
also of the solar zenith angle and the atmospheric conditions. 



ATMOSPHERIC EFFECTS ON RADIOMETRY FROM ZENITH OF A PLANE 
WITH DARK VERTICAL PROTRUSIONS 

1. Introduction 

In previous studies atmospheric effects on radiometric imaging from Landsat were analyzed. 
The effects are found to be quite small when three conditions are met: (i) the surface is 
bright, surface reflectivity larger than 0.20. (ii) the imaged area consists of large homogeneous 
fields, and (iii) the atmospheric absorption is small. This conclusion can be reached through 
radiative transfer calculations assuming an appreciable scattering optical thickness, as analyzed 
by Otterman and Fraser 1 . .Alternatively, the effects can be assessed resorting to a simplified 
radiative transfer treatment for an optically thin atmosphere, as formulated by Otterman et alv 
Lambert law was assumed for the surface reflection in both papers. This assumption is quite 
generally made in such studies. 

Another approach to this problem area is experimental in nature: a direct comparison of 
the satellite-measured surface-atmosphere system reflectivities with those of the surface ob- 
tained from field measurements. The comparisons involving different parts of the sandy 
Sinai Desert showed a close correspondence between the surface spectral reflectivities, deter- 
mined from the Landsat Multispectral Scanner System (MSS) digital data, with only a very small 
atmospheric correction, and those determined in the field using a hand-held Exotech -100 
radiometer 1,3 . The exception was the O.'-l.l nm MSS band, where a significant correction for 
water vapor absorption appeared appropriate. 

The Landsat measurements support the assumption of the Lambert law for the surface re- 
flection. inasmuch as the surface reflectivities of these sandy areas, bare of vegetation (northern 
Sinai and also Mexico near the Colorado River) derived at various times of the year (various 
solar zenith angles) are essentially constant: the variations in the reflectivity measured on dif- 

ferent dates arc small, and no systematic dependence in these variations with the solar zenith 
angle was reported 1 . 
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The task of monitoring a complex surface such as shown in Fig. 1 necessitates a quantative as- 
sessment of the structure of the clumps/plants that protrudes above the soil-plane, as well as 
a measurement of the soil-plane reflectivity. 

Our analysis here of this difficult radiometry is limited to conditions of a low optical 
thickness and of a laterally homogeneous atmosphere. Favorable conditions for remote sensing 
in and regions exist about 250 days a year. When the conditions are unfavorable, the aerosol 
optical thickness can be quite high. This, by itself, does not necessarily preclude classification 
of large fields by remote sensing from satellites 6 . However, associated with a large optical 
thickness, pronounced inhomogeneities are often encountered in the horizontal distribution of 
the aerosols. Such inhomogeneities can increase the difficulty in the quantitative remote sens- 
ing of the surface from satellites. 

We present first the previously developed surface model 5,7 aimed at describing the reflec- 
tivities of the terrain type shown in Figure 1. A treatment of the atmospheric effects over 
such a plane with sparse protrusions is then developed. Subsequently, these effects are analyzed 
in a comparison with those over a L. mbert plane. 

2. The Surface Reflectivities of a Plane with Sparse Protrusions 

The plants are modeled as thin vertical cylinders rising from the horizontal plane of the 
soil. The model is specified by only one geometrical parameter, s. which is the product of the 
height, diameter and number of cylinders per unit area (s is dimensionless). The soil and the 
cylinder mantles (the vertical area of the protrusions) are assumed as of a definite Lambertian 
reflectivity. r t for the soil interstices and r. tor the cylinders. The shadow cast on the hori- 
zontal plane of the soil by a single cylinder or several cylinders with a projection s on a verti- 
cal plane, is s tan 0 Q . It means that a fraction of 1-s tan 6 0 of a parallel beam at a zenith 
angle d Q is intercepted by the soil and a fraction s tan d 0 by the cylinders. This holds true 
as long as the shadow of one cylinder does not fall in part on another cylinder With lengthen- 
ing shadows, this assumption and the model become inaccurate. Eventually, at large solar zenith 
angles, the model breaks down. 



Overlapping of shadows can be taken into account if these fractional direct (luxes are ex- 
pressed by exponentials in -s tan $ Q , rather than hy the linear terms. The fractional direct 
flux is then exp (-s tan $ Q ) on the soil and 1 - exp(-s tan d Q ) on the cylinders. These expres- 
sions can be derived assuming a random distribution of cylinder-to-neighbor distances and azi- 
muths. This assumption should not be regarded as rigorously correct for any ecosystem, and 
especially not under the conditions of extreme aridity, where the requirement for a certain 
minimal plant-to-neighbor distance produces a fairly regular spacing of the protrusions. We 
stress thus the approximate nature of our model. In this study, we make one more approxi- 
mating assumption, that the protrusions are non-reflecting, that is, r c = 0. 

The surface reflectivity to the zenith accrues only from the interstices, inasmuch as the 
thin protrusions cannot be observed viewing to the nadir. In this direction the bidirectional reflec- 
tivity Tp for the direct solar beam or any beam at zenith angle d Q , is the product of the soil 
reflectivity ^ and the fractional flux on the soil plane: 

r p (V s ’ r i ) = r i exP^o) U) 

where rj 0 = tan d Q . 

Throughout the paper, the bidirectional reflectivity is formulated as ir times the radiance, 
resulting from an irradiance of 1 on a horizontal plane. In the case of a Lambert plane, the 
bidirectional reflectivity of the surface thus defined equals the hemispheric reflectivity (spectral 
albedo). There is no implication, certainly, that this holds for our surface. Thus. r p does 
not equal the hemispheric reflectivity, unless s * 0. The hemispheric reflectivity for a plane 
with sparse protrusions was derived and discussed elsewhere 7 . 

In field measurements of reflectivity, the reflectivity for the global irradiance is measured. 
The reflectivity r for the direct beam equals the field measurements only when these are con- 
ducted at a specific solar elevation, as discussed later. 
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3. The Reflectivity to the Zenith of the Surface-Atmosphere System for a Plane with Sparse 
Protrusions. 

3.1 Approach to the Calculation 

A plane parallel scattering atmosphere is assumed of a low vertical optical thickness r. r <<. 1 . 
The distribution with height of the scatterers is arbitrary. The simplified single scattering approx- 
imation (SSS) is used, as introduced into the same problem by Otterman and Fraser 8 : the num- 
ber of the scattered photons is exactly represented, but the second and higher order scattering events 
are disregarded. !n other words, the direction of a photon after the first scattering is frozen, and can- 
not be changed by subsequent encounters. This assumption applies to the photons scattered from 
the direct solar beam, as well as to those scattered from the flux reflected from the soil plane. 

A point-symmetrical phase function is assumed for the scatterers. that is, we assume an equal 
probability of scattering into each of any two opposite directions. The expressions obtained in 
the analysis apply thus directly to the Rayleigh as well as to the isotropic scattering. The down- 
ward directed pencils of the scattered radiation are represented by an equivalent flux. The magni- 
tude of the flux at the surface (above the protrusions) is represented as half of the total scattered 
from the beam, that is, as 0.5 [ 1 - exp (- t/m 0 ) 1 • where p Q = cos d Q .. The accuracy of this repre- 
sentation, analyzed elsewhere 9 , is satisfactory for our purposes. The average zenith angle of the 
flux scattered downward from the direct solar beam, calculated so that the illumination of the 
soil-plane is the same as by the scattered pencils, is separately analyzed 10 . 

We compute the radiance to the satellite at the zenith from an area within the field of view 
of its radiometer. We refer to this area as the object pixel. We initially assume that the surface 
is laterally uniform. In a later part we modify the solution to derive the adjacency effects of cross 
radiance and of cross irradiance. for a small object pixel embedded in a different, uniform, infinite 
terrain. 

The reflectivity to the zenith r n of a surface-atmosphere system can be expressed in a way 
that applies to a Lambert plane as well as to our complex surface, as 



r„ <0 o .r.s. r,) - r, 
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Gj* * 0 . r. s)(exp (-r) + F* (r. $)) 
_____ 


+ E(^ 0 . r), (2) 


where E expresses the baokscattering by the atmosphere to the satellite of the solar beam. Gj* is 
the fractional irradiance on the soil-plane (the interstices) by the solar beam and by the flux scat- 
tered from it. 2B** is the fraction of the flux reflected from the soil-plane that is backscattered 
by the atmosphere to the soil-plane. F* expresses the scattering to the zenith over the object pixel 
of the flux reflected from the entire soil-plane, and where exp<-r) represents the extinction of the 
radiance from the object pixel on the vertical path to the satellite. In this paragraph and through- 
out the paper, the word soil-plane is not synonymous with the surface: in ail the starred parameters, 
the probabilities of non-interception by the protrusions for the solar beam or for the pencils 
of the scattered or of the reflected radiation are incorporated. This non-interception has to be 
assessed 0 *: the downleg ir. G ( *. on the upleg in F*. and both on the upleg and the downleg 
in B**. Obviously, if s * 0. all the probabilities of non-interception are I and our equation 
reverts then to one for a Lambert plane with a reflectivity r . 

3.2 The Atmospheric Baokscattering Term. E(0 O . r). or Airtight 

The atmospheric baokscattering term E describes the radiance scattered to the satellite 
from the direct solar beam. This term is completely independent of the surface boundary, 
and can be referred to as the veil term. The solar irradiance on a horizontal plane at the top 
of the atmosphere is 1. The total of the (luxes scattered from the direct beam is I - e\p(-r ) 
where /a 0 * cos 0 O . Applying the SSS approach, the radiance to the zenith is expressed as 
( I - exp (-r u 0 ll Ptl80° - 0 O >. where the phase function p is normalized so that its integral over 
4t is 1. Multiplying this radiance by tr. as is our convention for the bidirectional reflectivity . the 
atmospheric backscattering or the veil term in Equation [2) is 


E t <? 0 . t ) - t| I -e\p(-r/i 0 >| Pi 180° — <# 0 » 


(3) 


o 
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For the limiting case (r!n 0 ) -*■ 0. E(0 O . r) * ir r P( 180° - 0 o )/m o - This expression becomes 
r/4/a 0 for the isotropic scattering and (3r/T6) 1(1 /m 0 ) + ,« 0 1 for the Rayleigh scattering. 

3.3 The Irradiance Gj*^, r. s) on the Soil-Plane. 

The irradiance by the direct beam above the surface (above the protrusions) is exp (- t/m 0 )» 
and the irradiance by the direct beam on the soil-plane is then exp(-r/M 0 > exp(-s tj 0 ). The 
downward-directed flux scattered from the direct beam is given by the SSS approach in the case 
of a phase function with point-symmetry as [ 1 - exp (~rin 0 )]!2. The fractional error of 
this representation of the downward flux is 0.15 r 2 !n 0 . 9 Only a certain fraction of this flux 
on the surface penetrates through the protrusions to the soil-plane. We denote this fraction 
e.xp(-$ rj x ). where rj x * tan 0 X . and $ x is the effective (average) zenith angle tor this penetrating 
fraction. When the distribution with the zenith angle of the irradiance L(0) cos 9 on a hori- 
zontal plane scattered from the direct beam is specified, the angle d x can be computed by 
integration: 

exp(-s tj y ) L(0) cos d sin 0d0 * *' * L(0)cos d sin 0exp(-stan0) d0 . (4) 

x 0 0 

The effective zenith angle 0 X depends thus not only on and the phase function of the 
scatterers. but also on the protrusion parameter s. For a phase function such as that 
in the Rayleigh scattering. 0 X changes only within narrow limits as a function of 0 o . 

We assume that 0, is independent of d Q . Variable absorption, in a layer close to the surface, 
such as by water vupor. can produce a marked change in 0 X . 10 In the numerical examples, 
below, we consider only one value of 0 X . such that rj x * 1.2. as applicable to a low value of 
s. in cases of a predominantly isotropic scattering and of significant absorption near the ground. 

An inappropriate selection of tj x causes only small errors when s is small. The problem vanishes 
if s * 0. 

The expression for the fractional solar irradiance on the soil-plane stands now as 
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G* ( 6 0 , r, s) » exp(-s tj 0 ) exp(-r/M 0 ) + 0.5 exp(-s tj x ) [1 - exp(-r/ju 0 )l 

« exp(-s n 0 ) |exp(-r/M 0 ) + 0.5 exp(-s(rj x - r/ 0 )) [I - exp(-r/M 0 ))} (5) 

Equation (5) in its last version is formulated to demonstrate that atmospheric scattering redirects 
a fraction 0.5 (1 - exp(-r/M 0 )J of the solar beam from a direction $ Q to the direction 6 % . In 
analyzing the atmospheric effects the consequences of this redirection are assessed. This redi- 
rection obviously does not affect the reflectivity in the case of a Lambert plane. The other 
half of the flux scattered from the direct beam is backscattered out of the surface-atmosphere 
system, thus reducing the irradiance Gj*. This reduction does not depend on the surface char- 
acteristics. 

3.4 The Fraction Backscattered to the Soil-Plane of the Flux Reflected from it. 

2B** (r, s), (Cross Irradiance). 

The factor 2B** describes the fraction backscattered to the soil-plane of the flux reflected 
from it. This factor carries a multiplier when contributing to the increase in the global ir- 
radiance Gj* of the soil plane, which is stated as 

G t * ($ Q , r, s) * Gj* ( 6 q , r, s) | 1 + 2r ; B**(r. s) + [2r> B** (r. s)] 2 . . . . J 

G*(d 0 , t, s) 

= ( 6 ) 

1 - 2rj B**(t, s) 

This simple formulation will not generally hold for an arbitrary complex surface: the flux 

backscattered to the surface will generally be distributed differently with the zenith angle 
(that is, effectively from a different direction) after each additional backscattering. and thus will 
have a different probability of penetrating to a specified surface element. Equation (6) applies 
for our plane with sparse protrusions, because the reflecting element is a horizontal Lambert plane 
and because the protrusions are assumed in our model to be located at random with respect to the 
illuminated/shadowed areas. Thus, the effective direction of a reflection does not depend on the 
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direction of illumination, and will be the same for any illumination of the soil-plane. 

The factor 2B** represents the probability of the atmospheric backscattering and also the prob- 
ability of non-interception by the protrusions on the upleg, for the flux reflected from the soil- 
plane, as well as for the downleg, for the backscattered flux. We introduce a factor B* such that 

exp (-s i ? r ) B* (r, s) « B** (r, s) (7) 

where 0 r • tan" 1 rj r represents the effective zenith angle for the interception by the protrusions 
of the backscattered flux. We assume later in the numerical examples 0 f * 0 X . 

In formulating the expression for 2B*, we consider a pencil of radiation cos 0/ir reflected 
from the soil-plane at a zenith angle 0. The zenith angle is denoted 0. 0 < 0 < ir/2. on the 
upleg as weil as on the downleg. The solid angle for all such pencils at an angle 0 is 2irsin0d0. 
The probability for these pencils to avoid interception by the protrusions and reach the atmos- 
phere is exp(-$tan 0). A nonintercepted pencil has a probability [1 - exp (— r/ cos0)]/2 of 
being backscattered to the surface (where the penetration to the soil-plane for the entire back- 
scattered flux is exp(-s rj r ). Summing the backscattering from all the nonintercepted pencils, 2B* 
can be stated as 

2B* (r. s) * cos0sin0 exp f-stan0) [1 - exp(-r/cos0)l d0 (8) 

0 

For the limiting case (r/cos0) -* 0 and for a Lambei. plane (s * 0), B* is the backscattering 
optical thickness r/2. As a numerical example consider r * 0. 1 , and a bright Lambert surface, r j * 0.4. 
The term 2B** constitutes then about 4% of the reflectivity r } . It is even a smaller fraction for a 
low reflectivity Lambert surface. For a plane with protrusions, the term ^ ZB** is reduced, due 
to the factors exp(-s rj r ), see Equation (7), and exp(-stan0) in the integral in Equation (8). 

3.5 Factor F*(r, s) for the Radiance Scattered to the Zenith Out of the Flux Reflected 
from the Soil-Plane (Cross Radiance) 

In Equation (2). ^ Gj*/(1 - 2r 4 B**) is the flux reflected from the soil-plane. In this 
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subsection the radiance scattered to the zenith over the object pixel out of this reflected flux 
is calculated, summing up the contributions from the reflection over the entire soil-plane. How- 
ever. we actually calculate the scattering to the zenith over the entire plane f rom the flux re- 
flected from a unit area. Then we involve the reciprocity law, to obtain the radiance to the 

zenith over the object pixel accruing by scattering from the pencils of radiation reflected from 
the entire plane. 

Consider a pencil of radiation r { cosd/tr reflected from a unit area of the soil-plane at a 

zenith angle 0. The solid angle for all the pencils at the same zenith angle is 2* sin0d0. A 

fraction exp(-stan0) of these reflected pencils will avoid interception by the protrusions. A 
fraction 1 - exp(-r/cos0) of a pencil that penetrates through the protrusions will be scattered 
by the atmosphere. The radiance to zenith will be a product of the irradiance of Gj/( l - 2rj B**). 
of alt the terms mentioned above and of P(0), the phase function for the angle between a pencil 
and the zenith. Therefore the parameter F* which describes this radiance in Equation (2). where 
it carries a multiplier r f Gj/fl - 2^ B**). is 

jt/2 

F*(t. s) * 2ir / sin0cos0exp(-stan0) P(0) [1 - exp(-r/cos0)l d0, (9) 

0 

where we have multiplied by ir. as is our convention in expressing the reflectivity r n in terms of 
the component radiances, see Equation (2). 

The radiance that we discuss in this subsection accrues from the reflection over the entire 
plane and affects the signal to the satellite over the object pixel. The phenomenon has been 
termed cross radiance by Otterman and Fraser 8 . It is further discussed in the su (section on 
the adjacency effects. We note that the protrusions (through the factor exp<-stan0) in the 
integrand) reduce the value of F* as compared to the cross radiance for the Lambert plane. 

The reduction is more pronounced for isotropic scattering than for the Rayleigh phase function, 
which peaks for 0 close to zero. In the abser.'e of protrusions, s * 0. and in the limiting case 
(r/cos0) -*• 0. F*(r. 0) * tv' 2 for cny point-symmetrical phase function. F*<0.1. s) is tabulated 
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for the isotropic and the Rayleigh phase functions in Table 1, for six values of s ranging up to 
0.30. For a limiting thin atmosphere, F*(r -*• 0. s) * rf(s) * B*(r -*■ 0, s) * rb*(s). 

3 6 Assessment of the Atmospheric Effects for a Uniform Plane 

We analyze the effects of a scattering atmosphere on the radiometric imaging of the sur- 
face, by comparing effects over a plane with sparse, dark protrusions, to those over a Lambert 
plane. It should be quite clear from our discussion that our results do not constitute a com- 
plete solution for complex surfaces where a significant reflection occurs from the vertical walls 
of the protrusions. Even though such reflected pencils do not directly reach a satellite at the 
zenith, atmospheric scattering from these pencils of radiation does affect the radiance received 
at the satellite. We first reformulate Equation (2) for the radiance to the zenith r n in terms 
of the reflectivity r p of our plane with protrusions, rather than in terms of the reflectivity r^ 

r n^o’ T> s> V ~ 

|exp(-r/M 0 ) + 0.5 f d (T? 0 , s) [1 - exp(-r/p 0 )]| (exp(-r) + F*(r, s)l 
r p (T? 0 , s, rj) ] _ 2r p (t 7 0 , s, Tj) f r (p 0 , s) B*(r. s) 

+ E(0 o ,r) (10) 

Equation (10) is identical to Equation (2), (in which Gj* is given by Equation (5) and B** is 
given by Equation (7), and where we substituted 

? d (rj 0 , s) = exp [-s(p x - r? 0 )l , (11) 

and 

f r (r? 0 , s) * exp [-s(77 r - r? 0 )] (12) 

We call f d and J r the redirecting factors. In what follows, we make a simplifying assumption 
that i? r = p x , and therefore the two redirecting factors are identical. f d ~ f r = f- For sun 
near the zenith, rj 0 < tj x , the redirecting factors are smaller than 1. and therefore the redirec- 
tion causes a reduction ir measured radiance r n . Conversely, when sun is at a low elevation. 
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t? 0 > i? x , the redirecting factors are larger than one, causing an increase in the soil-plane irrad- 
iance and therefore in the measured radiance (reflectivity). It follows from this discussion that 
the reflectivity rp for he direct beam is the same as for the global irradiance only when 
<9 0 * d x . An alternative treatment of the atmospheric effects, in terms of the reflectivity r g measured 
in the field at any solar zenith angle, is given in Appendix A. 

The differences r n - r p are plotted vs. r p for t 35 0.1 in Fig. 2 for s * 0.2 and s = 0.0. 

We refer to the r n - r p plots as Cs (complex surface) in the cases s = 0.2 and as La (Lambert 
plane) in the cases s * 0.0. An alternative abscissa, that of r i( is also indicated by marking 
each of the graphs at the points where r 4 = 0.1(0. 1)0.9. In the Figure, three pairs of plots 
are presented: Fig. 2A for tj q ■ 1.2 (d Q * 50.2°). that is, at the solar elevation at which the 

redirection factor f is 1 ; Figure 2B for t? 0 = 0.4 (ff Q = 21.8°), at which the redirection factor 
is smaller than 1, f * 0.85 if s = 0.2; and Fig. 2C for tj q = 2.0 (6 = 63.4°), at which the re- 
direction factor is larger than 1, f = 1.17 if s * 0.2. 

We first examine the plots of r n - r p for r? 0 = 1.2. f = 1.0. For r p = 0. the difference 
r n - r p stems from the atmospheric backscattering only (term E). and is therefore the same in 
the La and in the Cs graph. With increasing surface reflectivity, the difference r n - r p accrues 
(i) from the term E, (ii) from the extinction, the contribution of which to r n - r p can be re- 
garded as linearly proportional to the surface reflectivity, and (iii) from the multiple reflections 
from the surface, caused by the atmospheric backscattering; see Otterman et al. 2 for a fuller 
discussion of the Lambert-surface case. We concentrate here on the differences between r n - r p 
in the case of our complex surface and in the Lambert case for the same value of r p . The 
magnitude of this difference we denote as A(Cs, La). 

For small values of r p , the differences between the two graphs increase approximately 
linearly with r p . but remain small: A(Cs, La) < 0.0 1 for r p < 0.4. With increasing r p . the 
Cs graph descends more steeply from its E value (for r p = 0) than the La graph. This is due 
to the smaller value of F*(r. s) for s of 0.2 than that for a Lambert plane, s = 0. F* can be 
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regarded here as a partial compensation for the extinction on the upleg. In this context, the 
extinction on the path object pixel to satellite is described as 1 - exp(-r) - F* *> r - F*. The 
differences between Cs and La are here almost negligible if the comparison of the two graphs 
is made for the same values of r^ For r p higher than 0.4. the La graph curves upwards, re- 
maining above the -0.01 ordinate, whereas in the Cs graph the upward curvature is small, the 
graph dipping to the ordinate of about -0.03 for r p * 0.7 (rj = 0.9). The differences A(Cs, La) 
become larger than 0.02 in this region of r p . The large differences primarily stem from the 
fact that in the case of Cs plot, the 2fB* factor is very small, and therefore the upward curving 
of the plot is reduced as compared to the La plot. 

Case A. 6 q - 50.2°, has been selected as the first numerical example because the redirec- 
tion factor is 1. The reflectivity for the scattered irradiance then equals that for the 
direct beam, and therefore field measurements of reflectivity for the global irradiance determine 
r p . Only in those conditions, the difference r n - r p represents the difference between the 
satellite measurements and the field measurements. Case B. 0 o - 21.8°, has been selected be- 
cause tj q - tan0 Q * 0.4 is much smaller than in case A, and the redirection factor is significantly 
smaller than 1. f = 0.85. We note first in the Figure 2B that the value of E (that is, r n - r p 
for r p = 0) is only slightly reduced as compared to A. While ainmass for scattering from the 
direct beam is smaller by a factor of 0.69. the value of the Rayleigh phase function is signifi- 
cantly larger at the scattering angle of 158° as compared to that at 130°. (In the case of 
isotropic scattering, the value of r n - r p at r p = 0 would have been considerably smaller.) 

With increasing r p . the difference A(Cs, La) increases approximately linearly with r p , somewhat 
faster than in the case A. This faster increase is due to the smaller value of the redirecting 
factor, in reducing the enumerator of the expression for the Cs curve. The difference between 
the curves Cs and La is not very large for small values of r p : A(Cs. La) < 0.015 for r p < 0.4. 
For r p > 0.4 the gap between the two graphs opens quite wide. This is. as in case A. due to 
the upwards curving of the graph La and only insignificant curving of the graph Cs. In case B. 

the still weaker curving of the graph Cs is caused by the lower value of the redirecting factor 
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in the denominator of the expression for r R . 

In case C, 6 Q • 63.4°, rj Q = 2.0. the two graphs virtually coincide, that is, A(Cs, La) is 
practically insignificant. The similiarity of the atmospheric effects applies only if a comparison is 
for the same value of r p . This is due to the large value of the redirecting factor, f = 1.17. 
However, in the Cs graph, the values of r p are much smaller than the corresponding rj. Obviously, 
there are large differences between the two graphs if a comparison is made for the same value of r-. 

4.0 Adjacency Effects for a Small Pixel Surrounded by a Uniform Terrain 

In this section the atmospheric adjacency effects are analyzed for a specific situation, 
when a small area (the object pixel) is surrounded by a different terrain extending as a uniform 
area to infinity. Both the object pixel and the surrounding terrain are represented as a plane 
with sparse protrusions. By introducing s = 0, either one can become a Lambert surface. We 
recapitulate now our formulation of the surface-atmosphere system reflectivity, as consisting of 
four terms: veil, signal, cross irradiance and cross radiance. 

The veil term E(0 O , r), see Equation (3). is entirely independent of the surface character- 
istics. Thus, it can be defined as the reflectivity of the surface-atmosphere system if the sur- 
face is black (non-reflecting). The signal term constitutes the additional reading (that is. in 
addition to E) register ng at the radiometer if the object pixel does reflect (surface parameters 
r^ s), whereas the surrounding terrain is black. It is given as 

SR = r p (7? 0 , s, rj)|exp(-T/M 0 ) + 0.5 f d 0? o , s) [1 - exp(-r/M 0 )l J exp(-r) . (13) 

where r p is given by Equation (1). In this expression for the signal term it is assumed that 
the dimensions of the object pixel are vanishingly small as compared to the effective height of 
the scattering layer; the cross radiance and cross irradiance among the elements within the 
object pixel are neglected. 

The cross radiance term is the additional reading at the satellite (in addition to E) if the object 
pixel is biack, whereas the entire surrounding plane does reflect, characterized by and s. The 
cross raciance term is given as 
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r D (V *• |exp(-r//i ) + 0.5? d (r? , s) [1 - exp(-r/M 0 )l)F*(r. s) 

CR * y I L—l (14) 

1 - ? p (77 0 , s, Fj) f r (T? 0 , s) B*(t, s) 

where the barred parameters refer to the entire plane, which is assumed uniform, except for 
the small object pixel This explicitly presented term with F* is identical to the cross radiance 
CR as implicitly given (without the bars) in Equation (10) for s * s and Fj = r. 

The cross irradiance term results from a reflection by the surface within the object pixel 
of the additional irradiance, by the atmospheric backscattering of the flux reflected from the 


entire plane. The cross irradiance term is 

|exp(-r/M 0 ) + 0.5f d (r? o .s)( 1 - exp(-r/M 0 )l} r p (77 0 .s.? i )f r (n o ,s)B*(T^)exp(-r) 
1 - r p (i 7 0 , s. s) B*(r, s) 

(15) 


Cl = r (r? 0 . s, rj) ? r (r? 0 . s) { 


This complicated expression for the cross radiance is also implicitly included m Equation (10), 
where s = s and = r 4 . 

In what follows, the term r p f r B* in the denominator of both CR and CR is neglected. 

This reduces the values of CR and CR by only about 1%. (The fractional error in the differ- 
ence CR - CR can be larger, however). The simplified expressions for CR and CR are denoted 
as CR and CR. Our concern is an assessment of the fractional uncertainty or error in mea- 
suring the surface reflectivity, introduced by the adjacency effects. We refer to the expression 
for this uncertainty as the fractional difference in cross radiance. AFCR. It is obtained dividing 
CR - CR by the signal component SR. (see Equation 13): 


AFCR (0 O . t. s, s. ?j, Tj) 

? p (t? 0 , s, F|) jexp(-r//i 0 ) + 0.5f d (M o , s) [ 1 - exp(-r//i 0 )l 
r p<V s ’ r i ) i ex P(-r/A 0 ) + 0-5s* d (M o . s) (1 - exp(-r/M 0 )l 


= [C(d,., t, s. s. Tj, rp F*(r, s) - F*(r, s)] exp r. 


|- F*(r. s) - F*(r, s)| 


expr 


(16) 
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where the parameter 


C(0 O , r. s, s, rj, r t ) 


,(t? 0 , s, r^lexpC-T/Mp) + 0 5 f d (T? 0 , s) (1 « exp(-r/M 0 )l 


r p (i 7 0 , s. r^ |exp(-r/ju 0 ) + 0.5 f d 0? o , s) [1 - exp(-r/M 0 )l 


(17) 


is the effective contrast ratio between the surroundings and the object pixel. 

The factors controlling the adjacency effects are tabulated for various values of s in 
Table 1. These are: the cross radiance factors for a finite optical thickness for Rayleigh 
scattering 


F * (r, s) = — / [1 - exp(- — )1 exp(-s ^ 

m 2 0 fi M 


I? 3 3, 

— )( T + -M 2 )dM (18) 
4 4 


where the subscript m stands for molecular scattering, and for isotropic scattering. 


1 1 r \J\ - - 

B* (r. s) = — / [1 - exp(- ~ )J exp(-s ) d(i 

2 o M M 


(19) 


which is identical to the backscattering factor of the flux reflected from the soil-plane. Those 
two factors are computed for r* 0.1. We also present the factors 


such that for r 


f m '(s) - rf ' 0 


1 1 

b*(s) = — L 


exp(- s 
exp(-s 



3 3 , 

) (7 + 7^' ) d M 

4 4 



( 20 ) 

( 21 ) 


F m* < r > s >-" Tf m (s) 
and 

B* (r. s) r b*ls) 

Given the values of the cross radiance factors in Table 1 . we can readily assess the adja- 
cency effects in some numerical examples. In the case s = s, that is. when the protrusion 
parameter within the object pixel is the same as that of the surrounding terrain, the contrast 
ratio C is simply r i /r i . a contrast independent of s. d Q , or r. The expression for AFCR is 
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reduced to (C - 1) F* expr * (Fj - rj) F* expr/rj inasmuch as F* * F*. Consider a case Fj * 
0.3 and r t * 0.1, that is, a contrast ratio of 3, and assume Rayleigh scattering and r * 0.1. 

For s * s = 0.2. the adjacency effect AFCR amounts to 6.8%. Under the same conditions, 
but for a Lambert surface, AFCR amounts to 9.5%. We note that the adjacency effect AFCR 
is significantly reduced for a plane with sparse protrusions as compared to that for a Lambert 
surface. The reduction does not depend on the contrast ratio. For the value 0.2 of the pro- 
trusion parameter, the reduction is in the ratio F m * (0.1, 0.2)/F m * (0.1, 0.0) * 0.031/0.043 * 
0.72 or by 28%. 

The general case, when s ^ s, is more complicated, because the contrast ratio C as defined 
is a function of s, s, 6 0 and r as well as of Fj and r { . Under similar assumptions (Fj * 0.3, 

Tj = 0.1. Rayleigh scattering and r * 0.1). consider the case s » 0.2 and s * 0.0. The contrast 
ratio is 2.65 at d Q of 21.8°. 2.36 at 50.2° and 2.05 at 63.4°. The adjacency effect AFCR at 
these zenith angles is 4.2%. 3.3%, and 2.2%, respectively. At very large zenith angles AFCR 
reverses its sign. The results are quite different when a Lambert surface surrounds a plane 
with protrusions. If the reflectivities of the soil plane remain respectively Fj * 0.3 and r ; * 0.1. 
but when s * 0.0 and s * 0.2, the contrast C is 3 41 at zenith angle of 20.3°. 3.81 at 50.2° 
and 4.39 at 63.4°. The adjacency effect AFCR increases with the zenith angle, from 12.8% at 
20.3°, to 14.7% at 50.2° and to 17.5% at 63.4°. 

5.0 Discussion and Conclusions^ 

In this study we addressed the problem of atmospheric interference in radiometric imaging 
of Earth’s land surfaces from satellites. It is our assertion that in most cases land areas should 
not be represented as a Lambert surface and require a representation as a complex surface. A 
rather primitive model has been presented here: a Lambertian soil-plane from which dark 
(absorbing) protrusions arise. The dark protrusions do not contribute to the reflection from 
the surface. Their effect can be visualized as that of an absorbing layer of thin vertical needles 
or chaff that float above a horizontal Lambertian plane. Radiative transfer through such a layer 
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is treated by defining its optical thickness s as a projection on a vertical plane of needles over 
a unit area. The transmission through the layer as a function of the zenith angle 0 is given as 
exp(- s tan0). The extinction is zero on a vertical path and tends to infinity on a near-horizon* 
tal path. The effects as a function of 0 are thus significantly different than those of a plane- 
parallel layer of gaseous or aerosol absorbers. 

The radiative transfer through the scattering atmosphere is treated approximately, applying 
the simplified single scattering approach. The number of the scattered photons is computed 
exactly far the direct solar beam and for the radiation pencils reflected from the surface, un- 
der the assumption of a plane-parallel atmosphere. The approximation is only that a scattered 
photon retains the direction after its first scattering even if it undergoes subsequent scatterings. 

The effects of the atmosphere are here expressed as the difference between the reflectivity 
r n of the surface-atmosphere system measured by a satellite-borne radiometer at the zenith and 
the reflectivity to the zenith r p of the complex surface for the direct solar beam. We 
first report our findings for imaging large, laterally homogeneous areas. The differences r n - r p 
are found quite similar to the differences r n - a Q when imaging a Lambert surface with a 
reflectivity a Q , if the comparison is made for r p = a Q . This applies to low and moderate values 
of r p , r p < 0.4. For large values of r p . the difference r n - r p becomes more negative than 
the corresponding r n - a Q . This occurs because for a 0 >0.4, the plots r n - a Q vs a Q curve upwards, 
due to a term with a Q 2 that describes the second reflection from the surface (of the flux backscattered 
by the atmosphere). This effect, and thus the upward curvature of the r n - r p plots, is 
sharply reduced in the case of plane with dark protrusions. 

We compared the effects over a plane with protrusions, r n - r p . with those over a Lambert 
plane, r n - a 0 , on the basis of r p * a Q . Because of the protrusions, the reflectivity of the 
surface is lower than that of the soil-plane (the interstices). The reduction increases with an 
increasing solar zenith angle. Thus, we can state that the atmospheric effects over plane with 
protrusions with a reflectivity r ; of the interstices are really comparable with those over a 
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Lambert surface with a low reflectivity, approximately given as r { exp(- stan d Q ). For a low re* 
flectivity surface, the backscattering from the atmosphere to the satellite is a more important 
component of the r n - a 0 difference than for a bright surface. The backscattering to the surface 
of the flux reflected from it is of little importance. A plane with protrusions can be regarded in 
this context as a low reflectivity surface, with the reflectivity effectively decreasing with increas- 
ing solar zenith angle. Thus, our conclusions apply most forcefully when the sun is low, (see 
Figure 2C, where r p for s * 0.2 is much lower than r^). 

In imaging areas with dimensions not large compared to the effective height of the scat- 
tering layer, adjacency effects have to be considered. We analyze such effects only when the 
object pixel is a small area different from the surrounding terrain, that stretches to infinity as 
a uniform surface. This limiting case was analyzed previously for Lambert surface, retaining 
in the analysis only linear terms in the optical thickness 2 . For a plane with protrusions, the 
simplest case is when the protrusion parameter s within the object pixel is equal to that in 
the surrounding terrain, s. The adjacency effects in this case are significantly lower than in 
the case of a Lambert surface with the same surroundings-to-object-pixel contrast ratio and 
the same scattering optical thickness. However, when s # s, the situation is more complex. 

The magnitude of the adjacency effects can vary pronouncedly with the solar zenith angle, and 
even the sign can change, because the effective contrast surroundings-to-object-pixel depends on 
the direction of the illumination of the surface. Limitations of this analysis of the adjacency 
effects should be noted. When the cross radiance is computed as accruing from fluxes reflected 
up to a certain distance (rather that over the entire plane), the effect on the radiometry of a 
small field depends crucially on the shape of the phase function in the forward region. 11 

This study certainly is not the definite treatment of the atmospheric effects on radio- 
metric imaging of land areas. The subject deserves continued studies, applying the approach 
of modeling and Monte Carlo techniques, in addition to the programs of field measurements of 
the reflectivity. 
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Figure Captions 

Figure 1. A 1976 photograph of an exclosure in the northern Sinai, an area fenced off in 
1974. 

Figure 2. Calculated differences r n - r p for a plane with sparse protrusions, s * 0.2 (dashed 
lines), and for a Lambert surface, s * 0.0 (solid lines), at three zenith angles. 
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Figure 2. Calculated differences r n - r for a plane with sparse protrusions, s * 0.2 
(dashed lines), and for a Lambert surface, s = 0.0 (solid lines), at three 
zenith angles. 
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APPENDIX A 

The formulation of the atmospheric effects in terms of the reflectivity r p for the direct beam 
(that is, as r n - r p differences) has the advantage that r p is given by a simple explicit expression 
r p * Tj exp (-s tan0 o ), solely in terms of the surface parameters r 4 and s. It is this reflectivity r p for 
a direct beam that is measured in the laboratory, whether of an individual plant or an array of 
plants. The reflectivity measured in the field, which we denote r g , does not equal r p except when 
d 0 * 0 X . (From a practical viewpoint. r p and r g , are close enough in value to be considered equal if 
45° <0 o < 60° and the scattering optical thickness is not too large.) Thus, a formulation of the 
atmospheric effects as the r n - r g differences has the distinct advantage of comparing two measured 
quantities, one from the satellite and one from the ground. The significant disadvantage of this 
approach is that the expression for r g is more cumbersome than chct for r p . inasmuch as we formu- 
late r g as accruing from three separate streams (rather than .1 single stream in the case of r p ). These 
streams are: (i) the direct beam, m 0 exp (-r/*x 0 ). that penetrates to the soil-plane at 0 O ; (ii) the 
stream scattered from the direct beam. 0.5 n 0 [1 - exp(-r//i 0 )I , that penetrates effectively at 0 X ; 
and (iii) the stream Gj, scattered back to the surface from the flux reflected from it: 

2r* B*(r s) 

G b tfo.r.s.ri) = P 0 {exp(-ST? 0 )exp(-r//i 0 )+exp(-sr? x ) 0.5 [ 1 exp(-r/M 0 )l} T) (A 

that penetrates to the soil-plane at d r The expressions for the magnitudes of the two scattered 
streams are based on the simplified single scattering tSSS) approach. 

As a sum of the responses to these three streams, the surface reflectivity r g is a function of the 


atmospheric conditions as well as of the surface parameters and the solar zenith angle: 


r g (0 o , r, s, fj) 
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exp(-sn 0 ) exp(-r/M 0 ) + 0 5 exp(-ST? x ) [ 1-0.5 exp(-r/M 0 )l + exp(-sq r ) G b (fl 0 ,r,s, Tj) 
r * G O (0 O * r i> 

where the global irradiance G 0 on the surface at the top of the protrusions is the sum of the three 
streams discussed above. 

With these explicit relations between r g and the parameters describing the surface and the 
atmosphere, and applying the relations presented in the text between r n and the same parameters, 
the r n - r g differences can be studied. 


a-: 


(A-2) 


